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SUMMARY 

We have discovered that the isothermal equation of state for solids In 
compression Is of a simple, universal form. This single form accurately des- 
cribes the pressure and bulk modulus as a function of volume for Ionic, metal- 
lic, covalent, and rare gas solids. 


THEORETICAL ASPECTS 

The various classes of solids - covalent, metallic, Ionic, and rare gas - 
can be distinguished by their rather different types of Interatomic Inter- 
actions. These differences would appear to preclude a common description of 
the energetics of these systems. 

.' We have found, however, that equations of state (P, V, T, relations) for 
these classes of solids In compression are of a universal nature. In the fol- 
lowing we will first discuss reasons why one might expect to find a universal 
form for the equation of state (EOS) of solids. Next we will propose a uni-, 
versal expression. Finally, the considerable EOS data available will be used 
to test our proposed universal form. 

EOS- are valuable not only for prediction of thermodynamical properties 
but also to obtain Insight Into the nature of solid state and molecular 
theories (refs. 1 and 2). In reference 3 a general zero temperature EOS was 
reported which describes well the compression of metals (In the absence of 
phase transitions). This was derived from a cohesive energy which was found 
essentially to depend exponentially on the Interatomic spacing. That Is, It 
depends on the deviation of the Wlgner-Seltz radius from the equilibrium 
value, ry 5 - r^s e . (Here r w $ Is defined In terms of the volume V per atom 
as V = 4irr$s/3.) 
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Wave-function-overlap Interactions such as those occurring In metals are 
also found In varying degrees In covalent, Ionic, and rare gas solids under 
high pressures. There are additional types of Interactions Important In 
determining the EOS of Ionic and rare gas solids, however. Attractive cou- 
lomblc Interactions and attractive fluctuating dipole Interactions are also 
present. In addition. If we now consider nonzero temperatures, there Is a 
contribution to the pressure due to the variation of entropy. It Is of the 
form T(8S/3V)j for the Isothermal EOS. 

The pressure contribution due to electron overlap varies essentially 
exponentially with ryjs - r^s e . All other contributions to the pressure for 
solids on compression vary much less rapidly as a function of r^j. lor 
example, the entropy pressure Is T(3S/dV)y = aBT, where a Is the thermal 
expansion coefficient and B Is the Isothermal bulk modulus. The product <*B 
Is known (ref. 4) to be essentially Independent of V for T > 6 q, where eg 
Is the Debye temperature. Thus, the entropy contribution to the pressure Is 
expected to be nearly Independent of r^s, at least for 1 > 6 q. Also the 
Ionic or fluctuating dipole Interactions depend only on powers of r^s- Thus, 
these Ionic or fluctuating dipole Interactions provide a pressure component 
which Is also relatively slowly varying with r^s compared to electron over 
lap Interactions. 

Consequently, one may expect the shape of the EOS In compression to be 
determined essentially by the exponential wave function overlap. Note this is 
not true In expansion. For example, as the Interatomic separation becomes 
significantly larger than the equilibrium separation for Ionic solids, the 
wave-function overlap becomes small and the electrostatic Interaction between 
the Ions becomes dominant. This Is quite different from metals, where It has 
been demonstrated (ref. 3) that the exponential Interaction remains dominant 
well Into the attractive region. 


A NEW EQUATION OF STATE FOR SOLIOS 


■ From the foregoing discussion, we propose that an exponential expression 
similar to that found (ref. 3) for metals applies to all classes of solids in 
compression over a range of temperatures excluding phase transitions. This 
expression has the form 

P(X,T) = 3B t [(1 - X)/X 2 ]exp[n T (l - X)] (1) 


where B-] Is the Isothermal bulk modulus at zero pressure, X = (V/V q) 1/3 
(V 0 = 4irr$s e /3 Is the equilibrium atomic volume), and 

n T = 3[(3B/3P) 1 - l]/2 (2) 


Equation (2) follows from equation (1), where (3B/3P)y Is evaluated at zero 
pressure. To test this, we first rearrange terms: 

H ( X , 1 ) = ln[ PX 2 /3( 1 - X)] * lnB-j + n-,0 - X) (3) 

If our expression were universal, then one would expect Isothermal plots of 
experimental H(X,T) to be linear In 1 - X, with Intercept and slope yielding 
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By and (3B/3P)y, respectively. The effects of entropy and forces other, than 
those due to wave-function overlap enter only through these two parameters and 
v o- 


We have tested equation (3) against high pressure data from all classes 
of solids, excluding materials which have phase transitions within the experi- 
mental pressure range. We found It to be remarkably successful. In figure 1 
we have examples from each class of solid at a variety of temperatures: solid 

hydrogen (4.2 K data, ref. 5), solid xenon (150 K, ref. 6 ), cesium (295 K, 
ref. 7), molybdenum (293 K, ref. 8 ), sodium chloride (298 K, ref. 9, see also 
ref. 10), and magnesium oxide (293 K, ref. 8 ). 

Despite the fact that all classes of solids are represented from a vari- 
ety of experimental methods, the data are linear In (1 - X) as predicted by 
equation (3). This graphically Illustrates that EOS are, In fact, universal 
In form, our primary result. 

One may also find this way of plotting the data convenient for analyzing 
It. The values of B-[ and ( dB/aP ) 7 obtained from Intercept and slope of 
these plots as listed In the caption of figure 1 are In very good agreement 
with experimental values available In the literature (see, e.g., refs. 1, 5, 
and 11 ). 

One need only know two zero pressure equilibrium numbers, By and 
( 3B/3P ) 1 , at each temperature to obtain the pressure-volume curve through 
equation (3). Many of these equilibrium values are available In the litera- 
ture. For example, one can use ultrasonic values of By and ( 3B/3P ) 7 to 
predict P( V) through equation ( 1 ) and compare It with Independent high pres- 
sure experimental results. This Is shown In figure 2 for NaCl. NaCl was cho- 
sen because of the large amount of experimental data available which allows It 
to be used as a high pressure reference scale (refs. 1 and 10). Values of 
By = 2. 34x10^0 Pa and ( 3B/3P) y = 5.35 were taken from ultrasonic experiments 
(ref. 1). The agreement between the predicted pressures (eq. (1), solid 
curve), and the experimental data Is very good over a broad range of V/Vq. 

Another gauge of the accuracy of equation (1) Is found from a comparison 
with the Murnaghan EOS: 


P(X,T) = [B 1 /(3B/3P) t ] 


" -3(3B/3P) 
X 


1 


Equation (4) Is obtained by Integrating the linear approximation 


B(P,T) = B t + (3B/3P) t P 


(4) 


(5) 


so that equation (4) has increasing accuracy as P -» 0. The Murnaghan EOS 
Is perhaps the most extensively used EOS for solids and has been argued (refs. 

1 and 2) to be one of the "best" two-parameter EOS available. The dashed 
curve In figure 2 Is a plot of the Murnaghan EOS (eq. (4)). While good agree- 
ment with the data and equation (1) Is obtained at low pressures, the Murnaghan 
expression deviates from both at larger pressures. This Is a critical test of 
the universal expression equation ( 1 ), because It shows that the exponential 
form correctly describes the nonlinear pressure contributions which are not 
Included In equation (5). 
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A more severe test of nonlinear pressure effects Is shown In figure 3, 
where hydrogen data Is plotted for relative compressions of over 50 percent In 
volume. The data are due to Anderson and Swenson (ref. 5). The predicted 
pressures from the universal expression, (eq. (1)), are shown In figure 3(a) 
as the solid curve while the Murnaghan prediction, (eq. (4)), Is the dashed 
curve. The Input values used In both predictions are By = 1.7x10® Pa and 
( 3B/3P) y = 7.0, as listed In reference 5. Note these Input values are quite 
close to those obtained from the Intercept and slope, respectively, in 
figure 1. Again, the agreement between equation (1) and experiment is good, 
while at the higher pressures the Murnaghan prediction deviates signifi- 
cantly. Bulk modulus data are also given in reference 5, which provides 
another test. From equation (1) we have 

B(X,T) = (B t /X 2 )[2 + (n T - 1) X - n-^lexpU-jO - X) 1 (6) 

Figure 3(b) compares the experimental data (circles) with the predictions of 
equation (6) (solid curve) and of Murnaghan (eq. (5), dashed curve). Again, 
the prediction from the universal expression, (eq. (6)), agrees well with the 
data, while the linear (Murnaghan) prediction deviates at high pressures. 

We conclude that the EOS of all classes of solids in compression Is,, to a 
good approximation, a universal function. It can be represented by a simple 
exponential form which yields at each temperature accurate predictions of 
pressure and bulk modulus versus volume. 
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Figure 1. - Experimental values of H(X,T) = In CPX 2 /3(1 - X)] versus 1 - X, where x = (v/v 0 )^ and v Q is the equilibrium 
volume Der atom. The source of the data and values of B T and (dB /dP) T as obtained from the intercept and slope of the 
plots of figure 1 are respectively: hydrogen, ref. 5, 1.66x10° Pa, 7.33; xenon, ref. 6, 1.63x10^ Pa, 07; cesium, 
ref. 7, 1. 68x10^ Pa, 4.07; molybdenum, ref. 8, 2.66xlO n Pa, 3.99; sodium chloride, ref. 9, 2.35xl0 10 Pa, 5.15; 
and magnesium oxide, ref. 8, 1. 52X10 11 Pa, 4 81. 
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Figure 2. - Pressure versus reduced volume for NaCI. Solid curve, Eq.. (1). 
Dashed curve, Murnaghan expression, Eq. (4). 
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Figure 3. - Pressure (3a) and isothermal bulk modulus (3b) 
versus reduced volume for solid hydrogen. Experimental 
data is from ref. 5. Present model is Eqs. (1) and (6). 
Murnaghan expressions are Eqs. (4) and (5). 



1. Report No. 


2. Government Accession No. 


3. Recipient’s Catalog No. 


NASA TM- 87303 

4. Title and Subtitle 


5. Report Date 


Universality In the Compressive Behavior of Solids 


April 1986 

6. Performing Organization Code 


7. Author(s) 


506-53-12 

8. Performing Organization Report No. 


Pascal Vlnet, John Ferrante, John R. Smith, and 
James H. Rose 


9. Performing Organization Name and Address 

National Aeronautics and Space Administration 
Lewis Research Center 
Cleveland, Ohio 44135 

12. Sponsoring Agency Name and Address 


E-2702-1 

10. Work Unit No. 


11. Contract or Grant No. 


13. Type of Report and Period Covered 

Technical Memorandum 


National Aeronautics and Space Administration 
Washington, D.C. 20546 


14. Sponsoring Agency Code 


15. Supplementary Notes 


Pascal Vlnet, Visiting scientist from Laboratolre de Technologie des Surfaces, 
Ecole Centrale de Lyon, B.P. 163, 69131 Ecully Cedex, France; John Ferrante, 

NASA Lewi s\Research Center; John R. Sml th. General Motors Research Laboratories, 
Warren, Michigan 48090; James H. Rose, Ami's Laboratory - USD0E, Iowa State 
University, Ames, Iowa 50011. / \ | 



16. Abstract 


ll tm 


We have discovered that the Isothermar equation of state for solids in compres 
slon Is of a simple, universal form. This single form accurately describes the 
pressure and bulk modulus as a. function of volume for tonic, metallic, covalent, 
and rare gas solids. 


17. Key Words (Suggested by Author(s)) 

Equation of state 
Bulk modulus 
Compresslbl 1 ity 

18. Distribution Statement 

Unclassified - unlimited 
STAR Category 77 

19. Security Classif. (of this report) 

Unclassified 

20. Security Classif. (of this page) 

Unclassified 

21. No. of pages 

22. Price* 


Foi sale by the National Technical Information Service, Springfield. Virginia 22161 















National Aeronautics and second class mail 

Space Administration 

Lewis Research Center address correction requested 

Cleveland. Ohio 44135 

Official Business 
Penalty for Private Use $300 



Postage and Fees Paid 
National Aeronautics and 
Space Administration 
NASA-451 


IMASA 



